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PURPOSE:  The  Long-term  Fate  of  Dredged  Material  (LTFATE)  model  is  a  combined  local 
hydrodynamics  and  sediment  transport  model  used  to  determine  the  long-  and  short-term  stability 
of  dredged  material  mounds.  This  technical  note  (TN)  introduces  a  new  cohesive  sediment  transport 
submodel  for  LTFATE,  which  includes  a  combined  current-wave  shear  stress  calculation  and  a 
layered  sediment  bed  model.  LTFATE  can  be  accessed  via  the  Internet  at  the  following  address: 
http:/Avww.wes.amiy.mil/el/elmodels. 

BACKGROUND:  At  some  open-water  placement  sites,  transport  of  fine-grained  sediments 
outside  the  site  can  be  a  concern.  First,  fine-grained  sediments  often  are  cohesive  and  attract 
hydrophobic  contaminants  such  as  polychlorinated  biphenyls  (PCBs).  The  transport  from  dredged 
material  mounds  comprised  of  these  sediments  is  significant  because  of  the  potential  effects  of  the 
associated  contaminants.  Second,  even  uncontaminated  fine-grained  sediments  can  cause  negative 
biological  or  aesthetic  impacts  if  sufficient  quantities  are  transported  into  sandy  areas. 


Unlike  sands,  which  tend  to  settle  back  to  the  sediment  bed  rather  rapidly,  cohesive  sediments  remain 
in  the  water  column  for  significantly  longer  time  periods  and  therefore  can  be  transported  farther. 
This  can  result  in  significant  bed  contzunination  a  considerable  distance  from  the  original  mound 
location.  It  is  therefore  important  to  determine  the  stability  of  a  mound,  i.e.,  those  conditions  under 
which  sediments  in  a  mound  remain  in  place.  LTFATE  (Scheffner  et  al.  1995;  Scheffner  1996)  was 
developed  at  the  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES)  originally  for 
simulating  sand  transport  from  dredged  material  placement  sites.  Sand  often  is  used  to  cap 
contaminated  sediments.  However,  before  capping  occurs,  cohesive  mounds  are  exposed  to  the 
water  column  and  could  experience  sediment  transport.  In  addition,  due  to  the  additional  cost 
sometimes  associated  with  obtaining  sand  for  capping,  there  is  interest  in  using  suitable  fine-grained 
cohesive  sediments  as  cap  material  for  open  ocean  sites.  Therefore,  LTFATE  was  expanded  to 
include  fine-grained,  cohesive  silt  and  clay  transport.  However,  the  original  LTFATE  cohesive 
sediment  sub-model  was  basic  and  could  not  simulate  the  effects  of  large  storms  accurately.  Thus, 
the  improved  cohesive  sediment  transport  sub-model  described  in  this  TN  has  been  developed. 

Unlike  sands,  the  interparticle  forces  of  fine-grained  sediments  (due  to  their  small  mass)  are 
significant  when  estimating  transport  processes.  Sand  erosion  can  be  relatively  simply  related  to  the 
grain  size.  This  is  not  true  for  fine-grained  cohesive  sediments.  Several  factors,  including  grain  size 
distribution,  mineralogy,  bulk  density,  and  organic  content  have  been  demonstrated  to  significantly 
affect  the  erosion  rate.  In  fact,  sediments  that  at  first  glance  may  seem  similar  may  have  orders  of 
magnitude  difference  in  their  erosion  rates  (Lavelle  1984).  Because  of  the  various  processes 
influencing  erosion,  the  rates  tend  to  decrease  with  depth  below  the  sediment/water  interface  even 
for  sediments  of  consistent  grain  size  and  mineralogy. 
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This  TN  documents  the  layered  bed  processes  and  combined  current-wave  shear  stress  processes 
that  have  been  incorporated  into  the  cohesive  sediment  transport  sub-model  of  LTFATE.  Future 
TNs  will  address  additional  LTFATE  model  improvements,  including  bed  armoring,  combined 
sand/clay  sediment  bed  processes,  cohesive  sediment  flocculation,  and  cohesive  sediment  settling 
speeds. 

LAYERED  SEDIMENT  BED  MODEL 

As  previously  stated,  the  rate  and  method  by  which  cohesive  sediments  erode  depend  on  several 
factors,  including  grain-size  distribution,  orgamc  content,  pore  water  content,  and  mineralogy, 
among  others.  Therefore,  erosion  of  cohesive  sediments  varies  significantly  from  location  to 
location  as  well  as  with  depth.  A  commonly  used  method  to  relate  erosion  to  bottom  shear  stress 
has  been  incorporated  into  LTFATE.  This  method  relates  erosion  to  a  function  of  shear  stress  to 
some  exponential  power.  The  equation  for  the  erosion  rate  €  in  pounds  (mass)  per  square  foot  per 
second  is; 


A) 


(1) 


where  A)  and  m  are  site-specific  parameters,  (pounds  (force)  per  square  foot)  is  the  near-bottom 

shear  stress,  and  is  the  critical  shear  stress  below  which  no  erosion  occurs.  To  predict  erosion  at 
a  given  location  accurately,  several  values  of  Aq  and  must  be  used  in  the  vertical  direction. 
Therefore,  a  vertically  layered  sediment  bed  has  been  incorporated  into  LTFATE,  which  includes 
varying  values  of  Aq  and  x^r  These  parameters  can  vary  significantly  from  layer  to  layer,  and  erosion 
rate  experiments  on  sediments  extracted  from  the  site  are  the  best  method  for  determining  the 
vertically  varying  values.  If  such  data  are  not  available,  reasonable  values  from  similar  sediments 
should  be  used.  The  figure  below  provides  an  example  of  the  user-provided  distribution  of  these 
parameters  as  well  as  possible  initial  thicknesses  of  each  layer. 


Water  Column 

Sediment  Bed 

Xcf  =  5  X  10"*  Ibf/ft^ 

A)  =  8x10‘® 

Initial  thickness  =  0.025  ft 

Xcr=  1  X  10^  Ibf/ft^ 

Ao  =  4x10'® 

Initial  thickness  =  0.025  ft 

T«.=  5  X  10'®  Ibf/tt^ 

X 

II 

Initial  thickness  =  0.050  ft 

Tcf  =  1  X  10'^  Ibf/ft^ 

>lo  =  5x10'^ 

Initial  thickness  =  0.050  ft 

Tcf=2x10'^lbf/ft^ 

Ao  =  3.75x10'^ 

Initial  thickness  =  0.050  ft 

Xc^=  2  X 10'^  Ibf/ft^ 

A)  =  2.5x10'^ 

Initial  thickness  =  0.100  ft 

x«.=  2x10‘^lbf/ft^ 

A)  =  1  xio'^ 

Initial  thickness  =  0.100  ft 

Xcr  =  2  X  10'^  Ibf/ft^ 

Ao  =  1  xIO’^ 

Initial  thickness  =  0.100  ft 
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During  a  storm,  erosion  may  remove  any  layer  or  layers  completely,  leaving  only  the  bottom  layers 
exposed.  Initially,  the  top  layers  of  cohesive  sediment  have,  in  most  cases,  high  water  content  and 
are  easily  erodible.  Wave  and  current  action,  as  well  as  the  effects  of  benthic  organisms,  keep  these 
sediments  in  this  state  during  nonstorm  periods.  Bottom  sediments,  protected  from  the  wave  and 
current  action,  become  more  compact  and  are  more  difficult  to  erode.  Storm  action  will  remove  the 
upper  layers,  leaving  the  more  compacted,  erosion-resistant  sediments  exposed  to  the  water.  Also, 
freshly  deposited  sediments  are  added  to  the  topmost  layer  and  are  easily  resuspended.  Over  time, 
these  deposited  sediments  will  experience  natural  aging  and  compaction  by  freshly  deposited 
sediments  if  there  is  a  thick  enough  layer  to  protect  some  of  them  from  current  and  wave  agitation. 
The  compaction  and  aging  processes  are  not  incorporated  into  the  present  version  of  LTFATE.  The 
table  below  Qlustrates  the  erosion  rate  g  for  each  of  the  eight  layers  (if  that  layer  were  exposed  to 
the  sediment/water  interface)  for  conditions  of  a  13-ft  wave  in  82  ft  of  water  with  a  0.66  ft/s  steady 
current  in  the  same  direction  as  the  waves.  For  this  example,  the  value  of  exponent  m  was  set  to  2, 
which  experiments  have  demonstrated  is  a  reasonable  value.  It  can  be  seen  from  this  table  that  the 
erosion  rate  will  vary  by  orders  of  magnitude  between  the  top  and  bottom  layers. 


A) 
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5x10'^ 
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3.75x10'^ 

2  X 10'^ 

1.2x10"® 

2.5x10'^ 

2x10'^ 

7.8x10"^ 

1  X  10'^ 

2x10'^ 

3.1  X  10"^ 

1  X  10'^ 

3.1  X  10"^ 

As  previously  stated,  the  values  of  Aq,  m,  and  x„  are  considered  to  be  site-specific.  It  is  well-known 
that  the  erosion  rates  of  fine-grained  and  mixed  cohesive/noncohesive  sediments  are  related  to  the 
previously  mentioned  bulk  properties.  To  date  no  method  has  been  developed  to  associate  these 
bulk  properties  to  erosion  rates  despite  the  fact  that  these  properties  are  easily  measured,  i.e.,  no 
general  quantitative  theory  of  fine-grained  cohesive  and  mixed  cohesive/noncohesive  sediment 
erosion  and  resuspension  properties  exists.  Pending  research  under  the  DOER  Program  involves 
developing  methods  to  associate  erosion  rates  (and  vertical  variation  in  erosion  rates)  to  these  bulk 
properties  using  a  high-shear-stress  straight  flume.  The  methods  developed  by  this  research  will  be 
incorporated  into  LTFATE,  resulting  in  greatly  improved  ability  to  predict  fine-grained  and  mixed 
sediment  bed  erosion  during  storms  and  a  useful  tool  for  improved  site  management.  Prior  to 
availability  of  these  erosion  rate  methods  related  to  bulk  properties,  this  flume  is  being  used  to  obtain 
high  shear  stress  erosion  rates  for  specific  sites.  This  information  is  being  incorporated  into 
site-specific  applications  of  LTFATE.  Current  studies  include  the  Los  Angeles/Long  Beach  Harbor 
and  New  York  Mud  Dump  sites. 
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COMBINED  CURRENT/WAVE  SHEAR  STRESS 

The  method  incorporated  into  LTFATE  to  estimate  bottom  stresses  due  to  combined  current  and 
wave  action  (Xf^  in  Equation  1)  is  described  in  this  section.  Large  water  bodies  such  as  lakes,  oceans, 
and  estuaries  are  subject  to  both  current  velocities  and  wave-generated  orbital  velocities.  As  a  result, 
the  bottom  shear  stress  will  depend  on  the  relative  magnitude  of  both  velocities  and  the  angle 
between  them.  In  general,  the  bottom  shear  stress  is  assumed  to  follow  a  quadratic  law: 

P/c“' 

2  (2) 

where  x^b  is  the  near-bottom  shear  stress  due  to  currents  (Ibfi^ft^),  p  is  the  density  of  water, /<.  is  the 
current-related  friction  factor,  and  u  is  the  current  velocity  outside  the  boundary  layer.  For  pure 
currents, /c  ranges  from  0.002-0.005,  depending  upon  the  bed  roughness.  If  waves  are  present,  a 
similar  equation  for  the  wave-related  shear  stress  Xy^,|,  is  assumed: 

-t 

2  (3) 

where  U  is  the  near-bottom  orbital  velocity  for  the  wave  and/,„  is  the  wave-related  friction  factor. 
The  coefficient/,,,  ranges  from  0.002-0.05,  depending  on  U  and  the  wave  period 

Often,  for  simplicity,  these  two  effects  are  simply  added  together.  As  will  be  described  later,  this 
will  produce  a  shear  stress  that  is  too  low,  often  by  a  factor  of  two  or  more. 

The  physical  problem  is  described  by  Grant  and  Madsen  (1979).  Summarizing  the  main  points  of 
the  wave-current  interaction,  the  near-bottom  flow  is  mfluenced  by  the  relatively  low-fi:equency 
currents  and  high-frequency  surface  waves.  The  bottom  boundary  layer  is  assumed  to  consist  of  an 
oscillatory  wave  boundary  layer  nested  within  a  relatively  steady  current  boundary  layer.  When  the 
maYimiim  bottom  orbital  velocities  of  the  waves  are  of  the  same  order  as  the  steady  current 
velocities,  the  small  scale  of  the  wave  boundary  layer  causes  the  boundary  shear  stress  that  would 
be  associated  with  the  wave  to  be  much  greater  than  that  associated  with  the  current  alone.  The 
wave  and  current  interact  to  generate  a  shear  stress  that  is  different  from  that  generated  by  the  sum 
of  the  two  components.  This  interaction  is  a  nonlinear  process  and  numerous  assumptions  are 
necessary  to  derive  it 

The  process  discussed  here  relies  upon  the  same basic  physical  description  as  that  provided  by  Grant 
and  Madsen  (1979)  but  is  somewhat  simpler  in  its  description  of  bottom  shear  stress.  Bottom  shear 
stresses  are  predicted  using  a  method  developed  by  Chnstoffersen  and  Jonsson  (1985),  referred  to 
hereafter  as  CJ.  Shear  stresses  predicted  by  the  CJ  model  compare  well  with  experimental  data. 
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As  defined  in  the  CJ  model,  bed  shear  stress  due  to  combined  current  and  wave  action  may  be 
calculated  from: 


^  r  2 
/\P  J  w  ^wbm  ^ 

2  (4) 

where  is  the  maximum  bed  shear  stress,  and  u^bm,  the  amplitude  of  the  bottom  orbital  velocity 
at  the  top  of  the  wave  boundary  layer  (or  maximum  bottom  orbital  velocity),  is  calculated  from 
linear  wave  theory  by: 


4ncosh(hH,) 

where  is  the  wave  height,  g  is  the  acceleration  of  gravity,  k  is  the  wave  number,  T,  is  the  wave 
period,  and  A  is  the  total  water  depth. 


m  is  calculated  by: 


m  =  (l+a  ^  +2alcos(5  -a)|)”^ 


(6) 


'^cb  _ 

\fc] 

'  u  ' 

^wbm 

L/J 

_^wbm_ 

(7) 


where  a  is  the  ratio  of  x^b  (bottom  shear  stress  due  to  currents)  and  the  amplitude  of  (Le., 
x^bm  maximum  value  of  the  oscillatory  x^^,  bed  shear  stress  due  to  waves).  5  is  the  angle  of 
the  current  direction  and  a  is  the  angle  of  the  wave  direction.  Note  that  if  x^b  » '^wbm^ 


m— » 


\fc] 

’  U  ' 

UJ 

J^wbm^ 

(8) 


and  thus  Xf,^  Xcb-  Similarly,  if  x^b  « then  m  1  and  X£,„  x^bm- 

Ultimately,  the  prediction  of  the  shear  stress  depends  on  the  turbulence-related  friction  factors. 
There  arc  separate  friction  factors ^  and/^  for  current  and  wave-related  processes.  It  is  beyond  the 
scope  of  this  TN  to  describe  the  methods  for  estimating  the  friction  factors.  It  will  suffice  to  say  that 
the  friction  factor  equations  are  experimentally  determined.  The  procedure  for  computation  of 
bottom  shear  stress  is  as  follows: 

a.  Compute /c  for  pure  currents  by  the  following  equation  (note,  for  pure  current  the  second  term 
on  the  right-hand  side  of  this  equation  is  zero): 
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'2’ 

f 

'30h 

In 

-In 

Uj 

K 

K 

(9) 


where  kff  is  the  Nikuradse  roughness  coefficient  (assumed  to  be  0.015  ft),  is  the  apparent 
roughness  (described  below),  k  is  the  von  Karman  constant  (=0.4),  and  h  is  the  water  depth. 

b.  Compute and  parameter  J  for  pure  waves  (m=l): 


/.=2 


.  ^wbm  . 


(10) 


^wbm 

■tfl 

J^N^a  . 

L  2  J 

(11) 


where  is  the  maximum  value  (amplitude)  of  the  oscdlatory  wave  particle  velocity  at  the  top 
of  the  wave  boundary  layer,  ^  Oie  angular  frequency  and  P  is  an  experimentally 

determined  constant  (=0.0747).  Now  recompute/^,,  using: 

2Pm 

(12) 


c.  Keeping/c  fixed,  iterate  to  find  a,  m,  /,  and/^  using  Equations  6, 7, 1 1,  and  12. 

d.  Compute  5,^  (wave  boundary  layer  thickness),  and  kj^  (apparent  roughness): 


2 


|0^ 


(13) 


*;A  =  305^exp 


V"iJ 


(14) 


where  r  is  an  experimentally  determined  constant  (=0.45). 

e.  Recompute^  from  a.;  this  time  the  second  term  on  the  right-hand  side  will  not  be  zero. 

f.  Repeat  steps  b.  through  e.  until/^  converges. 

g.  Calculate  x^,  maximum  bed  shear  stress,  using  Equation  4. 
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Now  that  the  maYimiim  bottom  shear  stress  due  to  both  currents  and  waves  is  known,  erosion  rates 
can  be  estimated  as  described  in  Equation  1. 

WES  is  involved  in  several  projects  where  the  improved  LTFATE  model  is  being  applied.  These 
projects  include  erosion  rate  parameter  information  from  the  high  shear  stress  flume  for  better 
quantification  of  the  vertical  variation  of  erosion.  A  future  TN  will  provide  details  on  the  results  of 
these  projects.  These  applications  will  be  relevant  to  site  management  concerns. 

POINTS  OF  CONTACT:  For  additional  information  on  LTFATE,  contact  the  author  of  this  TN, 
Dr.  Joseph  Gailani  {60l-634-4S5lJ.gailani@cerc.wes.army.mil)  or  the  DOER  Program  Managers, 
Mr.  E.  Clark  McNair,  Jr.  (601-634-2070,  mcnairc@exl.wes.army.mil)  and  Dr.  Robert  M.  Engler 
(601-634-3624,  englerr@exl.wes.army.mil).  This  technical  note  should  be  cited  as  follows: 

U.S.  Army  Engineer  Waterways  Experiment  Station.  (1998).  “LTFATE  cohesive 
sediment  transport  model,”  Technical  Note  DOER-Nl,  Vicksburg,  MS. 
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